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We study the adsorption of ideal random lattice copolymers with correlations in the sequences on 
homogeneous substrates with two different methods: An analytical solution of the problem based on 
the constrained annealed approximation introduced by Morita in 1964 and the generating functional 
(GF) technique, and direct numerical simulations of lattice chains averaged over many realizations 
of random sequences. Both methods allow to calculate the free energy and different conformational 
characteristics of the adsorbed chain. The comparison of the results for random copolymers with 
different degree of correlations and different types of nonadsorbing monomers (neutral or repelling 
from the surface) shows not only qualitative but a very good quantitative agreement, especially in 
the cases of Bernoullian and quasi-alternating random sequences. 



I. INTRODUCTION 

Random copolymers (RC) attract much attention and have been studied intensively in the last 
two decades. This is motivated, on the one hand, by the nontrivial properties of individual RCs 
and their behaviour in solutions and melts, and, on the other hand, by the biological relevance of 
this class of molecules (protein [l[ and RNA Q folding, hybridisation of heterogeneous DNA 0, 
because real proteins and nucleic acids represent, with no doubt, irregular copolymers. 

Among the many open questions involving RCs, the problem of RC adsorption at different in- 
terfaces is taking a particularly prominent place, because of the technological relevance of such 
systems, their role as biomimetic model systems for molecular recognition Q etc.), and the ba- 
sic challenges they offer for theoretical physicists and applied mathematicians. Several types of 
interfaces have been considered, interfaces that are penetrable for polymer such as liquid/liquid 
interfaces or lipid membranes, and impenetrable solid substrates, either with or without chemical 
or geometrical homogeneities S S 0,11, I, M El IH El- In the present paper we focus on the 
problem of adsorption of a single ideal RC chain with correlations in the monomer sequence on an 
impenetrable planar surface using a lattice model. 

Solving problems that involve RCs, one faces the necessity to carry out two operations: Firstly, 
one needs to sum over the statistical weights of polymer conformations, i.e. to calculate the 
system's partition function Z and the free energy F (the logarithm of the partition function, 
F = —ksTlogZ), and secondly, one needs to average the free energy with respect to the disorder 
in the RC monomer sequence. Averaging "log Z" over the sequence disorder turns out to be 
a very difficult task which can seldom be carried out rigorously, ft would be much easier to 
average the partition function prior to taking its logarithm, but this corresponds to the physically 
different situation of annealed disorder (in our system the monomer sequence is quenched), which 
approximates the quenched situation only very poorly. One common approach consists in using 
the replica trick which then can be combined with some variational scheme (e.g., the reference 
system approach suggested by Chen @, H, H , or the Gaussian variational approach Q). 

In the present work, we exploit another approximation to resolve the difficulties with the 
quenched average - the constrained annealing suggested by T. Morita in 1964 14]. In this ap- 
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proximation scheme, one averages over annealed disorder (i.e. directly averages the partition 
function), but imposes the constraint that the first, the second, etc., moments of the monomer 
distribution in the RC sequence keep their correct values. The ideas of the method are developed 
in the original paper of Morita a more detailed and formalized description can be found in 
the paper of Kiihn [15| and in a recent review of Soteros and Whittington [16| . 

For the problem of RC adsorption on the solid interface and RC localization at the liquid-liquid 
interface, the Morita approximation has been successfully applied by Whittington and Orlandini et 
al. [111. 1 1 21] . In these works, the authors considered simplified polymer models - fully directed ran- 
dom walks in 2D space (Dyck and Motzkin paths). RC localization at penetrable and impenetrable 
interfaces [ll[ (with both chain ends fixed at the interface) as well as the adsorption-desorption 
transition upon pulling one chain end away from the interface by an external force [l2| were stud- 
ied. The Morita approximation was first order, i.e., only a constraint on the first moment of 
the monomer distribution was imposed. (Note that the annealed approximation can formally be 
considered as a zero order Morita approximation.) 

Higher order Morita approximations were applied to the RC adsorption problem by Alvarez et 
al. [13| . where polymer trajectories were again modeled by Dyck and Motzkin paths. The authors 
obtained lower bounds limiting the quenched free energy and demonstrated how these bounds 
improve as the order of the Morita approximation increases. In contrast, the location of the phase 
boundary is not more accurate in the Morita approximation than in the annealed approximation, 
as was rigorously proved by Caravenna and Giacomin [TtJ 

The simplified polymer models (Dyck and Motzkin paths) considered in the previous works [TTL 
H3 . HH have their doubtless advantages of being exactly solvable and at the same time permitting to 
obtain a physically clear and tractable picture of the studied phenomena. However, real polymers 
"live" in the 3-dimensional space and are not necessarily directed. Furthermore, the RCs in Refs.fTll. 

were of Bernoullian type, i.e., uncorrelated. In the present paper we study - although 
somewhat simplified - a 3-dimensional situation and a wider class of RCs with correlations in the 
monomer sequence, where Bernoullian RCs represent one special case. 

The calculation of the partition function via summation over polymer conformations is carried 
out using the generating function (GF) approach (or grand canonical approach) developed for 
polymeric systems by Lifson [l8l ]. This very general method has been used to study conforma- 
tional transitions in polypeptides (coil - a-helix and coil - /3-structure [2(|) and DNA (coil 
- double helix) [2ll |22| . Birshtein was the first to apply it to the problem of single polymer ad- 
sorption [23|, |24| . She developed a theory of the adsorption transition based on the GF formalism 
and demonstrated its power and generality with a number of different examples. The approach 
allowed to reproduce very elegantly a number of previously known results, and, in the same easy 
manner, to treat the more complicated case when adsorption is coupled with a helix-coil transi- 
tion in the chain. In her work Birshtein [23l |24| focused on the adsorption-desorption transition 
point and the order of the transition. Later, the results of her theory were used to analyze the 
adsorption-desorption transition for regular (A m B n ) x -multiblock copolymers J25[. Brun 12(1 has 
proposed an approximate generalization of the equation for the transition point in Ref. |23j for 
correlated annealed RCs (see section UlI A[) . 

The method of GF is highly general and its success is due to the possibility to calculate the 
GF for elements of the chain conformations. Despite the fact that it is more than 30 years old, 
it is still often used up to date. One of the most prominent examples is the ongoing discussion 
about the order of the DNA denaturation transition [2?], [H, [2i| [3(| , which is carried on - for the 
most part - in the framework of this approach. We also note that in its "continuous" version, 
the method of GF is equivalent to the propagator formalism in Fourier-Laplace space - see, for 
instance, Muthukumar et al. 

In the present work, analytical calculations based on the Morita approximation and the GF 
approach (though the final equations can only be solved numerically, the disorder average and the 
summation over chain conformations are carried out analytically) are supplemented with numerical 
simulations for tethered RC chains. In the simulations, a sample of several (sufficiently many) RC 
sequences are drawn randomly from the desired distribution. For each chain, the propagators are 
calculated exactly (recursively). This allows to evaluate the free energy and different observables 
directly and carry out an "exact" quenched average over the sequence disorder. The comparison 
of the results obtained with the two different methods allows us to make a conclusion about the 
accuracy and applicability of the Morita approximation. 
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FIG. 1: Random copolymer near adsorbing plane. 



The rest of the paper is organized as follows. Section |TT] defines the model and introduces the 
annealed and the Morita approximation. Although the annealed approximation is deficient in 
many respect, it is convenient to introduce the GF approach in this framework. This is done 
in Section III CI An additional motivation for treating the annealed case in detail is that there 
is a class of so called two-state or annealed copolymers, which are adequately described by the 
annealed approximation (see Ref. [33[ and references therein). The generalisation to the Morita 
approximation in Section III CI is quite simple and transparent. In Section III Dl our numerical 
approach is introduced. The results obtained using both methods are compared and discussed in 
Section ITTT1 Finally, Section HVl presents the summary and final conclusions. 



II. MODEL AND METHOD 



A. Definition of the model 



Consider a RC chain consisting of N monomer units near an impenetrable surface, Fig. [T] A 
lattice model of the polymer is employed, i.e. each chain conformation is represented as a walk 
on the lattice. Monomer units in the RC chain can either be of type A or of type B, and their 
sequence is taken to be random. We assume that there are no excluded- volume interactions between 
monomers (i.e., we consider the case of ideal chains), but there is a short-range interaction between 
monomers and the surface. The Hamiltonian of the system can be written as 

JV 

H = A * ' bti*A + (i - xM • (i) 

i=l 

Here \i indicates the type of the monomer i: 

{1 for A-monomer (giving [...] = £^4) , , 

for B-monomer (giving [...]= Eb) 

and Aj is a conformation-dependent parameter, namely, Aj = 1 if i-th monomer occupies a site 
adjacent to the adsorbing surface and otherwise. Any contact of an A monomer with the surface 
leads to an energy contribution ea, and likewise, a contact of a B monomer contributes Eb- 

The monomer sequence \ = {xiiX2, ■■■iXn} is random and will be represented as a first order 
Markov chain. It is determined by the probabilities to find A and B monomers in the sequence 

P{X 1 =A)=f A and P{ Xl = B) = f B = l-f A (3) 

and by the probabilities that the monomer of the type i is followed by the monomer of the type j, 

P (Xrn =j\Xm-l = = Pij 

(4) 



( PAA 


PAB \ 


1 = 1 


( PAA 1 - PAA 


V PBA 


PBB J 




V 1 - PBB PBB 
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The transition probabilities are normalized : . Pij = 1 . We also assume that the Markov chain is 
reversible, i.e. fiPij = fjPji- Hence, the probability of the sequence \ = {xi,X2, -,Xjv} to appear 
is given by the product 

P(X) = P(xuX2,---,Xn) = f Xl -P X1X2 -P X2X3 ■ ••• -Pxn-l,Xn- ( 5 ) 
It is convenient to introduce the cluster parameter c: 

c := I - pab - Pba = Paa+Pbb - l, (6) 

which characterizes the correlations in the sequence: c > means that there is a tendency in the 
sequence for grouping similar monomers into clusters, c < favors the alternating sequence of A's 
and B's, c = corresponds to uncorrelated (Bernoullian) sequences. The cases c — —1 is also 
special, since it describes perfectly ordered alternating block copolymers {...ABABAB...). 

Because of the normalization and reversibility conditions, it is sufficient to prescribe any two 
of these parameters to completely determine the first order Markov chain, for example, Ja and c. 
Then 

PAA = !A + C(l - f A ) , PAB = (1 - /a)(1 - c) 

Pba = /a(1 - c) , PBB = 1 - /a + cf A - 

Finally, it should be noted that the range of possible values for /a depends on the value of c: 
for c > Ja may take any value between and 1 whereas for c < there is a constraint 

<f A <^-- (8) 



1-c ' 1-c 
In the actual calculations of this paper, we consider two cases: 

• A monomer is attracted by the surface (sa = —1), B monomer is neutral (e_e = 0), denoted 
as SN (sticky-neutral) case 

• A monomer is attracted by the surface {ea = —1), B monomer is repelled from the surface 
(sb = 1), denoted as SR (sticky-repulsive) case 

It should be emphasized that the terms "neutral" or "repulsive" are applied to monomer-surface 
interaction whereas a polymer chain composed of either neutral or repulsive monomers will effec- 
tively repel from the surface. 



B. Disorder average and Morita approximation 

Having introduced the model, let us now write the system's partition function for a given real- 
ization of the monomer sequence 

Zn{P\x) =]Texp{-/3]>>. [ Xl e A + (1 - Xi)e B ]\ , (9) 

w I »=1 J 

where /3 is the inverse temperature and ^ w denotes the sum over all chain conformations oj. 

Rather than dealing with particular realizations of a RC sequence, we study the whole ensemble 
of RC chains with the desired statistics. The randomness in the RC sequence is quenched, which 
means that after synthesis of a RC chain, its monomer sequence remains unchanged. To obtain the 
quenched free energy, the logarithm of the partition function should be averaged over all possible 
sequence realizations: 

(3F q = - (log Z N ((3\ X )) (10) 
where the angular brackets (. . .} denote averaging over sequence randomness 



(A) :=Y,P(X)A( X ). 

x 



(11) 
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Direct averaging of the log Zjy is, in general, a rather difficult problem, so one can use, for example, 
the replica trick and work formally with an m times replicated system 0, H US EE H3] ■ 

Averaging of the partition function prior to taking the logarithm corresponds to the annealed 
type of disorder. In this case, the average is evaluated according to 

pF a = -\og{Z N ((3\ X )) 1 (12) 

where 

(Z N (0\x)}=J2 P (x)^expl-f3j2 A i-iXieA + (l-Xi)eB}Y (13) 

X uj { i=l ) 

Physically, annealed disorder means that the type of any monomer (and its affinity to the surface, 
in our case) may change while the system attains its equilibrium state [33[ . The annealed approx- 
imation can roughly be considered as the "zero order" approach to the problem of carrying out a 
quenched average. It is substantially easier to perform but it often fails to give good values for the 
free energy and other observables. The annealed approximation does not guarantee that even the 
lower moments of the distribution of monomers remain correct in the final result. 

To approach this problem, Morita [14j | has suggested to use a partial annealing procedure where 
the moments of the monomer distribution are constrained to keep their correct values. For example, 
in our case, the constraint on the first moment reads 



^5> = /a, ( 14 ) 

i=l 

the second moment of the distribution is constrained according to 

,Xi-iXi = fAA, (15) 



N 

V 



N-l^>- 

i=2 



where Jaa — j A • PAA is the AA dyad probability, and so on. For the averaging of the parti- 
tion function, this means that we introduce the constraints (fT4|) - (fT5|) into Eq. (fT3| via Lagrange 
multipliers A and k: 



(z n ((3\x))m = E p w E ex p | -0 E A * • ^ £A + ^) £s l f 

x exp | A (j2xi - Nf A ^j + K (^Xi-xXi - (N - l)f AA J j 



(16) 



and request 

d\og(Z N ((3\ x )) M 
d\ 

d\og(Z N ((3\ x )) M 



N 



3 k 



J2xi)-Nf A = 0, 

I N \ y J 

z2x<-iX>)-(N-l)fAA = Q. 



\i=2 



Here and below the notation (. . .) M for the disorder average subjected to the Morita conditions 
is used. After exchanging the order of conformational and disorder average (i.e. the order of 
summation J2 x J2lo = X^S X ); ^ ne partition function reduces to the sum of the pre-averaged 
statistical (Boltzmann) weights over polymer conformations. In the present work this sum is 
calculated using the GF approach (or grand canonical approach) developed for linear polymers by 
Lifson O. 
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C. Sum over conformation and GF approach 



The main object (concept) of the GF approach is the GF 

oo 

r(a:) = ]T Z N x N , (18) 

N=l 

where Zn is the partition function of a polymer chain with N units. The GF can be considered 
as the grand canonical partition function; in this case, x = exp(/i/fcsT), where /i is the chemical 
potential of a monomer unit and x is called monomer activity. Once the GF |T5J) is known, the 
partition function Z^ can be formally determined as the coefficients of the expansion of T(x) in 
powers of x. In the long chain limit, N 3> 1, the asymptotic expression for Zn is even simpler 
since it is dominated by the smallest singularity x c of T(x): Zn ~ x~ N . 

The GF approach is particularly efficient in the cases where monomer unit may exist in one 
of several states (for example, double helical or coil in DNA [111, [22||, a- helical and coil or 
/3-structure and coil [20| in polypeptides) and the chain conformation can be represented as an 
alternation (not necessary regular) of sequences of different types. 

To study RC adsorption, we generalize the work of Birshtein, who was the first to apply the GF 
formalism to the problem of homopolymer adsorption [23l l2~ij . For clarity, we first illustrate our 
approach for the case of annealed RCs. 

For an adsorbed chain, each adsorbed conformation can be represented as a sequence of adsorbed 
and desorbed groups of monomers. The adsorbed sequences are called trains, among the desorbed 
we distinguish between tails at the ends of the chain and loops separating adsorbed sequences. It 
is clear that an adsorbed chain can have at most two tails, i%t < 2, whereas the number of loops is 
less than the number of trains by unity, i.e. if nz = ns — 1- We note that different tails, loops, and 
trains do not interact with each other. The GF of the partition function, averaged over (annealed) 
sequence disorder, T(x) = 53jv=i (Zn(P\x)) is then given by 

T(x) = [FP^TTix) + (Wf)TR- 1 ] • T s (x) 

x [E-TL^Tsix)]- 1 • [T T (x)+E)-e. ' " 

Details of the calculation are given in the Appendix [XI Here f = (/a, fs) T is the vector of single 
monomer probabilities, P is the transition probability matrix, Eq. the matrix R is defined as 
R = PW with the "interaction matrix" 

W = ( W Q w B ) ' WA = ' WB = e ~ (iSB ' (20) 

E is the unity matrix, and e is the vector e = (1, 1) T . The functions Ts(x), Fl(x), and Tt{x) are 
the GFs of adsorbed sequences, loops, and tails, respectively, in matrix form, 

oo oo 

r s (a;) = «sW(PW)V = n S (n)R n x n 

n=l n=l , . 

oo oo V / 

T L (x) = n L (n)P n x n , and T T (x) = £ Q T (y)P n x n , 

n—l n—1 

where fi^ (n) is the number of conformations corresponding to each sequence of the type i of length 
n. We note that T(x) is still a scalar function. The expression for r(x) has the same structure 
as that obtained by Birshtein for homopolymer adsorption in Ref. [231 ] . except that it is a matrix 
equation in the RC case. 

To find the smallest singularity of T(x), one needs to consider the smallest singularity associated 
with [F, -T L (x)T s (x)}~ which can be found as the smallest root of the equation 

det[E -T L (x)T s {x)} = 0. (22) 
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It must then be compared with the smallest singularity xy of the GF for the free chain in a bulk 
(in the absence of a surface), 



IV (a:) 



oo 

E 



n v (n) x r ' 



(23) 



Since the GFs of matrix arguments in Eq. (f22|) are matrix power series, they can be easily calculated 
using eigenvalues and eigenvectors of the corresponding matrices [35| (P and R). Despite its 
simplicity, Eq. (|22"|) can be solved only numerically. The only case when it can be simplified 
corresponds to Bernoullian, i.e. uncorrelated, copolymer, where c = 0. Applying the above 
scheme, it is straightforward to show that the determinant equation (|22p reduces to the scalar one 



1 - T L (x)Ts(xw e ff) = 



(24) 



The effective w e jf — Jawa + (1 — /b)i«b, in accordance with earlier finding [lq . 

The method of GF allows the calculation of various properties of the adsorbed chain. One of 
the most important quantitative characteristics of polymer adsorption is the fraction of monomers 
in the surface layer. Using the definition of (Z n (j3\x)), Eq. ([9]) and differentiating x c with respect 
to wa and wb, one can calculate the fraction of A- and B-contacts with the surface: 



<IA.B 



d log x c 



WA.B 



dX c 



d((3e A ,B) x c dwA.B 

The derivatives of x c with respect to wa,b are obtained from the derivatives of lhs of Eq. (f2"2"|) 



(25) 



dXr 



8D/dwA,i 



dw A , 



dD/dx 



where we have defined 

D :=dct [B-T L {x)Ts(x)]. 
To differentiate the determinant D we use Jacobi's formula 

ddetX(r) {>,., v , dX 

— — = Trace Adi (X) • — — 

dr \ dr 



(26) 



(27) 



(28) 



where "Trace" stands for the trace of a matrix and Adj(X) denotes the adjugate matrix for X. 
Hence, 



3D 

dx 



= — Trace 



Adj(E -T L (x)T s (x)) 



dx 



r s (a;)+r L (a:) 



dT s 

dx 



and 



3D 

dw A ,i 



= — Trace 



Adi(E~T L (x)T s {x))-T L (x) 



dT s 

dw A ,_ 



(29) 



(30) 



After having derived the appropriate formalism for annealed RC adsorption, we are now ready 
to generalize our results to the Morita approximation. The latter differs from the annealed approx- 
imation in that it imposes conditions on the first and second moments of monomers' distribution, 
Eq. (fT"7|) . Comparing the expressions for the partition function, Eqs (fi"3"|) and (Tl6|) . it is easy to see 
that the calculation of the GF in the Morita approximation T(x) = {Zn{(3\x)) m xN ^ s ec l urv " 
alent to that in the annealed approximation, if the following changes are made: x — » xe 
and 



-A/a — k/aa 



VAAe X+K PAB 
PBAe X PBB 



(31) 
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The matrix P is again R = PW. If we denote y :— xe X ^ A K ? AA , we obtain the analog of Eq. 

det[B-T L (y)T s (y)}=0. (32) 

To find the values of the variational parameters A and k, the Morita conditions (JTTJ) are applied. 
In terms of y they can be rewritten as 

TTT h J A — U , 

J X (33) 
d}ogy 

— H 1- JAA — U 

OK 

where the logarithmic derivatives of y can be expressed by differentiating Eq. (f3"2")l with respect to 
A and k. By solving equations (j3"3"|) together with Eq. (|3"2"|) , the values of y c , A and k and, therefore, 
the smallest singularity x c = y c e x ' A+K ^ AA , as functions of (3, ea and are obtained. 

The calculation of observables is similar to that described for the annealed case with the only 
change x c — > y c . Indeed, x c = x c [wa, wb, M w a, wb), k(wa, wb)) and 

d\ogx c _ d\ogx c d log x c d\ d log x c dn 
dwA dwA dX dwA dn Owa 

d\ogx c , d\ogx c , . . , . 

Since — — — and — are equal to zero because of the Morita conditions, we get 

o\ an 

dhgXc = dlogxc = d\ogy c ^ 
dWA dwA dwA 

and the same is, of course, valid for the — — -. 

dwB 

The final missing ingredient in the theory is the actual form of the GFs of loops and adsorbed 
sequences, rg(x) and Tl(x). They depend on the details of the particular system under considera- 
tion, i.e., on the lattice type and on the geometry of the adsorbing substrate. Here we consider the 
case of the adsorption on a plane, when polymer conformations are simple random walks, restricted 
to the half-space, on the conventional 6-choice simple cubic lattice (SCL). 

Calculating the GF for adsorbed sequences is a relatively simple task. Using the definition of 
the GF |18]), we obtain: 

r s (a:) = x + (z - 2) x 2 + (z - 2f x 3 + . . . = X (36) 

1 — (z — 2)x 

Here z = 6 is the coordination number of the SCL. 

The calculation of the GF for loops is more complicated. It is described in Appendix [Bl The 
result is 

T L {x) = i- {l - (z - 2)x - y/(l - zx)[l - (z - 4)a;]} . (37) 



D. Numerical simulations 



Our analytical calculations are complemented by numerical simulations of RC adsorption on 
planar surfaces. The numerical method is described below. The model is the same (sectioning] but 
the system under consideration is slightly different. First, the analytical calculations were carried 
out for chains of infinite length, whereas in the numerical simulations, the chains must remain 
finite for obvious reasons. The main consequence is that we no longer have a sharp adsorption 
transition, but a smooth crossover region between an adsorbed and a desorbed regime. We have 
not studied chain-length effects in this work. Second, the chains in the numerical simulations are 
tethered to the surface at one end. This helps to avoid an uncertainty with the normalization, 
because for free chains of finite length one would have to introduce a box of finite size. We note 
that tethered chains and free chains have the same adsorption characteristics in the infinite-chain 
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limit, hence our analytical results also apply to tethered chains (in our theory, the main equation 
(|22p or Eq. only contains GFs for adsorbed sequences and loops and is therefore equally valid 
for both free and tethered chains). 

The numerical treatment is based on a Green's functions formalism first introduced by Rubin 
[3rl l37l . [HI and later used in the more general theory of Scheutjens and Fleer [3§| . The central 
quantities are the statistical weights Gt(z;n) of all conformations of tethered RC chain parts of 
length n with one free end in the layer z, and the corresponding weights Gf(z;N — n) of all 
conformations of chain parts between the nth monomer and the end monomer N, subject to the 
constraint that the monomer n is fixed in the layer z whereas the end monomer is free. They 
satisfy the recurrent relations 

G t j(z; n + l) = {\G tJ (z - 1; n) + (1 - 2X)G tJ (z; n) + XG t . f (z + l;n)} (38) 

at z ^ and 

G t (0; n + 1) = exp(-/3e n+ i) {(1 - 2A)G t (0; n) + XG t (1; n)} , (39) 



G/(0; n + 1) = cxp(-/3 £A ,_„ +1 ) {(1 - 2A)G f (0; n) + AG t (l; n)} (40) 

with e„ = Xn£A + (1 — Xn)^B, where A is the probability that a random walk step connects 
neighboring layers. On the simple cubic lattice, one has A = ^. Using as the starting point the 
monomer segment distribution G t (0;l) = exp(— /fei), G/(0;1) = exp(— an d G t (z;l) = 0, 
Gf(z; 1) = 1 for z > 0, and recursively applying Eqs. (f3"8T[4T)|) , one can easily calculate Gtj(z; n) for 
every chain length n. The statistical weight of all conformations of tethered chain is then obtained 
by summing over all positions of the free end: 

Nl 

Z{N)^Y. G ^ N ) ( 41 ) 

2=0 

The change in the free energy of the tethered chain with respect to the free chain in the solution is 
given by AF = — log Z{N). Here the translational entropy of the free chain has been disregarded, 
i.e., the chains are assumed to be sufficiently long that it can be neglected. At the transition point 
one has AF = 0, i.e., the energetic benefit of monomer-surface contacts is equal to the entropic 
penalty caused by tethering the chain to the plane and thereby restricting the number of possible 
conformations. 

Using Gt and Gf one can also calculate the probability that monomer n is adsorbed on the 
surface via the composition law 

This gives the average fraction of adsorbed monomers 

JV 



q N 

71=1 

and the average fractions of A and B contacts with the surface 

1 N 1 N 

qA = Nz2p( n ) S XiA > <lB = jj^2p(n)5 Xt , = q-q A . (44) 

n—l n—1 

In this paper, we present results for chains of length N = 1000. We consider two-letter correlated 
RC sequences with the distribution defined in Section III Al For every set of model parameters Ja 
and c, the adsorption characteristics were calculated as a function of the inverse temperature /3 
for 100 different sequence realizations and then averaged. This corresponds to a situation with 
quenched sequence disorder. 
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III. RESULTS 



A. Annealed chains 



We first discuss the adsorption properties of annealed chains, focussing on the results for the 
two cases: sticky-neutral (SN) and sticky-repulsive (SR) defined in Sec. Ill Al 

Fig. [2] shows the results of the numerical solution of Eq. (f22|) as a function of the inverse 
temperature (3 at Ja — 0.5 for various values of c in the SN and SR cases. One can see that 
blockier RCs with higher values of c have lower x c (hence, a lower free energy (3F/N = logx c ). 
At small values of [3 (at high temperatures), all the curves terminate in the adsorption-desorption 
transition point where x c = xy (recall that xy is the smallest singularity of Yy(x) - the GF for 
the bulk chain, Eq. (|23[) ; for a simple cubic lattice one has xy = 1/6). The transition point, the 
critical value (3t r of (3, is therefore found from the equation 

det[V -T L (xy)T s {xy)} = 0. (45) 

The resulting value of (3t r as a function of /a, the relative weight of A- monomers in the sequence, 
is shown in Fig. [3] For c < the allowed range of /a narrows according to Eq. ([8]) . In the case 
c = — 1 it degenerates to a single point Ja — 0.5. It can be seen that the aggregation of monomers 
into clusters shifts the transition points and the x c ((3) curve towards lower values of (3 (i.e., it 
favors adsorption) and this shift is more pronounced in the SR case. 

The equation for the transition point (|45[) has some similarity with the result obtained 
by Brun [261 ]. However, the latter contains only the largest eigenvalues, v m ax, Umax, of 
the matrices P and R, respectively. In our notation the equation of Brun has the form: 
1 — T L(v max xy)T s{HmaxXv) = 0. It gives the same result as Eq. (|45l) in the Bernoullian case 
(c = 0), but overestimates f3 tr at c < and underestimates it at c > 0, compared to our result. 

Figs d] and O show examples of the temperature dependence of the fraction of A- and B-contacts 
for the two cases SN and SR. 

At c = —1, where the chain is a regular alternating copolymer, the "annealed approximation" of 
course becomes exact and the results are in excellent agreement with the corresponding numerical 
results, which are also shown for comparison. The remaining difference is due to the finite chain 
length in the numerical calculations. This situation can be used to illustrate the main qualitative 
difference between the two cases SN and SR: The fraction of B-contacts (Figs. [4]b and[5]b) behaves 
differently depending on whether B-monomers are neutral or repelling. In the former case qs grows 
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FIG. 4: Fraction of A and B contacts in the annealed approximation for /a ~ 0.5, sa = — 1, £b — and 
different values of cluster parameters c = — 1, —0.75, —0.5, —0.25, 0, 0.25, 0.5, 0.75 as indicated (lines). 
Also shown are numerical results for c = — 1 (crosses). 



with (5 whereas in the latter case it behaves nonmonotonically. The initial growth is caused by 
a cooperativity effect: B-monomers linked to adsorbed A-monomers have a higher probability to 
come into contact with the surface (see Fig. [6]). In the SN case, such contacts are not penalized, 
hence the B-monomers are effectively attracted by the surface and qs continues to grow for all (3. 
In the SR case, the cooperativity effect is counterbalanced by the increasing effect of the energy 
penalty on B-contacts, hence qs decreases again at high (3. 

Next we consider the curves corresponding to truly disordered chains with c > — 1. Both in the 
SN and SR cases the fraction of A-contacts grows with increasing f3 up to qA = 1, which is much 
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FIG. 5: Fraction of A and B contacts in the annealed approximation for /a = 0.5, sa = — 1, £b = 1 and 
different values of cluster parameters c = —1, —0.75, —0.5, —0.25, 0, 0.25, 0.5, 0.75 as indicated (lines). 
Also shown are numerical results for c = — 1 (crosses). 
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FIG. 6: Illustration of cooperativity effect: Adsorbed train and its adjacent monomers. 



higher than the expected fraction j a of A-stickers in the sequence (not necessarily adsorbed) for 
quenched disorder ((?a(/3 3> 1) < /a=0.5). At the same time, the fraction of B-contacts exhibits a 
nonmonotonic behaviour both in the SN and the SR case: Close to the transition point it grows, 
and then, after passing through a maximum, it decreases and virtually vanishes. This reflects the 
physical nature of the annealed approximation mentioned above: Any monomer can change its 
affinity to the surface (from negative to positive), if it is thermodynamically preferable for the 
system, and this is exactly the case at high (3. Close to the transition point, the degree of A^B 
conversion is small and neutral or repelling B-monomers come into contact with the surface due to 
the cooperativity effect described above; of course, for given chain statistic (given f a and c) and 
given inverse temperature /3, qs is smaller in the SR case than in the SN case due to the additional 
penalty on B-contacts. 

This example demonstrates that the annealed approximation gives us essentially incorrect result 
regarding one of the most important conformational characteristics of adsorbed heteropolymer 
chains. Therefore it is necessary to improve on this approach. The Morita approximation (Sec. Ill 5)) 
provides us with such a possibility. 



B. 'Morita chains' and Quenched chains 



The objectives of the present work are not only to apply the Morita approximation for the 
problem of adsorption of correlated RCs, but also to estimate how good it can approximate the 
quenched (i.e. the real) system. The results obtained using the Morita approximation and the 
numerical approach (Sec. Ill D[) are thus presented together. We will use the following convention: 
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FIG. 7: Fraction of A and B contacts calculated in the Morita approximation (lines) and using the transfer 
matrix approach (symbols) for /a = 0.5, ea = — 1, £b = and different values of cluster parameters 
c = -0.75 (□), -0.5 (O), -0.25 (O), (A), 0.25 (v), 0.5 (+), 0.75 (★). 



The results obtained with the Morita approximation are shown as solid lines, whereas the data 
for the corresponding quenched systems are represented by symbols. Since our primary interest 
is to study the influence of correlations and the nature of the non-adsorbing monomer B on the 
RC adsorption properties, we set /a = 0.5 in this section. The results for c = — 1 are not shown 
because they coincide with those obtained using the annealed approximation and have already 
been presented in Figs 0] and [5l 



C. SN case 



We consider first the SN case. Fig.[7]shows the temperature dependence of the fraction of A and 
B contacts for different values of the cluster parameter c. Both qA and qs grow with increasing (3. 
The results for qs differ qualitatively from those obtained for all random copolymers (c < —1) in 
the annealed approximation, but they resemble those obtained for alternating copolymers (c = —1) 
in Fig. |U The reason for the monotonic growth of qs is the absence of an energetic penalty on a B 
contact, combined with the cooperativity effect discussed earlier (Fig. [5]). This effect explains the 
larger amount of B-contacts for RCs with lower c at intermediate and strong adsorption (far from 
the transition point). Contrary to the annealed case (Fig. |Ja), the value of qA does not exceed 
fx- 

Comparing the results obtained using the Morita approximation and the numerical procedure, 
we have not only qualitative, but also a quantitative agreement. The curves for qA are predicted 
almost perfectly by the Morita approximation, the results are only very slightly overestimated. 
The qB curves are in reasonable agreement for c < 0, whereas for c > 0, the Morita approximation 
underestimates qg. This can be explained by the fact that the Morita copolymer has the possibility 
to adjust itself to the surface - A and B monomers can be redistributed between trains and loops. 



D. SR case 



Let us turn to the SR case. The results for the monomer-surface contacts qA and qs are presented 
in Fig. [5] As in the SN case, the fraction of A contacts grows with increasing (3. Furthermore, qA 
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FIG. 8: Fraction of A and B contacts calculated in the Morita approximation (lines) and using the transfer 
matrix approach (symbols) for /a = 0.5, ea = — 1, es = 1 and different values of cluster parameters 
c = -0.75 (□), -0.5 (O), -0.25 (O), (A), 0.25 (v), 0.5 (+), 0.75 (★). 



decreases with decreasing cluster parameter c, in a much more pronounced fashion than for the 
SN case. In addition, the growth of qA with becomes slower with decreasing c. 

Unlike qA, the fraction of (unfavourable) B contacts qs depends non-monotonically on the in- 
verse temperature, Fig.[5]b. It increases just after the transition point, passes through a maximum 
and then decays. This behavior looks similar to that observed in the annealed approximation (both 
for the SN and SR cases, Figs. 2] b and[5]b), but its nature is very different: In the annealed case, 
B-monomers can convert to A ones, whereas both in the quenched case and in the Morita approxi- 
mation a B— >A conversion is not possible. The initial growth of qs is due to the cooperative effect, 
Fig. [6l whereas the subsequent decay is caused by the increasing influence of the energy penalty 
on B-contacts. (Indeed, one can easily imagine the most probable (ground state) conformation in 
the j3 — ► oo limit: All A-monomers are on the surface, all B-monomers are desorbed in loops.) It 
is worth noting that the Morita approximation seems to underestimate qA and overestimate qs 
at large /3, i.e., small temperatures. This is counterintuitive, given that the Morita copolymer 
can adjust the A-B distribution on the sequence to the conformation to some extent. In fact, the 
discrepancy is simply an effect of the finite chain length in the numerical simulations. On the one 
hand, the cooperativity effect is reduced in the free tail, which reduces the number of B-contacts. 
On the other hand, A-monomers close to the tethered end can make contact with the surface more 
easily than A-monomers in the middle of the chain, which enhances the number of A-contacts. 

The results presented in Fig. [5] lead us to the conclusion that the Morita approximation repro- 
duces the fractions of A and B contacts in the quenched system reasonably well. The quantitative 
agreement is best for quasi-alternating RCs (c = —0.5 and —0.75). We suggest the following ex- 
planation of this effect: The Morita approximation reduces the initial heteropolymer problem to 
a homopolymer problem, thus "smoothing out" the monomer sequence of the RC. Obviously, a 
quasi-alternating RC can easily (and more realistically) be represented as a homopolymer with a 
monomer unit replaced by an effective one than a RC with c > where the blocky structure is 
essential (it could be better represented by a multiblock copolymer). 



E. Adsorption transition point in the Morita approximation 



Solving the problem of RC adsorption in the Morita approximation amounts to finding the 
solution of a system of three equations, Eqs. ((32)) and (f33|) . with three unknowns: y c , A, and n (the 
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partition function/the smallest singularity of T(x) is then given by x c — y c e x ^ A+K ^ AA ). 

We obtained that the Lagrange multipliers A abd k vanish in the the adsorption transition 
point: \((3tr) — «(/3tr) — 0- This implies that at the transition point, the Morita and the annealed 
approximation are equivalent. Hence, the transition point is the same for the Morita copolymer 
than for the annealed copolymer and is determined by Eq. (I45|) . consistent with the conclusion of 
Caravenna and Giacomin 11711 . 



IV. SUMMARY AND OUTLOOK 



In the present work, the adsorption of a single two-letter RC chain with correlations in its 
monomer sequence onto a planar homogeneous surface was studied. One of the monomers (A) was 
always taken to be attracted to the surface whereas the other (B) - was neutral or repulsive (SN: 
sticky - neutral or SR: sticky - repulsive cases, respectively). The monomer sequence of the RC 
was modeled as a first order Markov chain, which is completely determined by two independent 
parameters. For the sake of convenience, we used as these parameters the probability for A 
monomers in the sequence, Ja, and the cluster parameter c (c = - uncorrelated (Bernoullian) , 
c > - quasi-blocky, and c < - quasi-alternating copolymer) . To average the free energy over the 
sequence disorder , i.e., over the whole ensemble of monomer sequences, two approximations were 
employed: The annealed approximation and the Morita approximation. The partition function 
was calculated in the grand canonical ensemble using the GF approach which is based on the 
representation of an instantaneous conformation of the adsorbed chain as a sequence of independent 
non-interacting trains (adsorbed sequences) and loops flanked by tails. 

In the framework of the annealed approximation, where the partition function, instead of its 
logarithm, is averaged over sequence disorder, a general approach for the calculation of the av- 
eraged partition function of the RC chain (in the form of the smallest singularity of the GF), 
the fraction of adsorbed monomers, and the transition point was developed. It was found that 
the annealed approximation cannot describe the quenched system correctly in the intermediate 
and strong adsorption regimes: For example, it strongly overestimates the fraction of contacts of 
adsorbing monomers, because of the possibility of A— >B monomer conversion. Nevertheless, the 
results for annealed chains are interesting in their own right since they describe the adsorption 
behavior of so-called two-state polymers (annealed copolymers). 

The theory developed for the annealed approximation was then extended and combined with the 
Morita approximation. The latter approach is based on a constrained annealing procedure, where 
the quenched average is replaced by an annealed average subject to constraint on the moments 
of the monomer distribution. In our case, constraints on the first and the second moments were 
imposed. 

To estimate the accuracy and validity of the Morita approximation, numerical calculations using 
a transfer matrix approach were carried out for the quenched system: Polymer chains of length 
N = 1000 monomers were considered, for the given monomer statistics 100 realizations of random 
sequence were generated, for each realization the free energy and other observables were calculated 
and then averaged over the different realizations. The comparison of the results obtained using 
the Morita approximation and those for quenched system showed not only a qualitative but also a 
good quantitative correspondence, especially for quasi-alternating RCs (i.e., RCs with c < 0). 

It was shown that the fraction of adsorbed A-monomers grows monotonically as the inverse tem- 
perature (3 increases, independently of the character of the interaction between the nonadsorbing 
B-monomer and the surface. In contrast, the fraction of adsorbed B-monomer behaves differently 
in SN and SR situations: In the former case, it increases with /3, whereas in the latter case, 
it behaves nonmonotonically, exhibiting first a relatively rapid growth just above the adsorption 
transition point which is followed by decay at higher values of (3. 

The results of our calculations also indicate that both the annealed and the Morita approxima- 
tions give the same adsorption transition point. 

The special case = 0.5, c = — 1 corresponds to the limiting situation of a quenched regular 
alternating AB-copolymer which can be solved exactly using the formal framework of the annealed 
approximation. 

We remark that the approach developed in the present work can be applied in a straightforward 
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way to other types of lattice models of polymers or to other geometries of adsorbing substrates. 
The key difference will be the calculation of GFs for adsorbed sequences (which is usually quite 
easy) and loops (which is a more sophisticated task). In the framework of the GF approach other 
interesting observables can be calculated as well, for example, the statistics of trains and loops. 
This will be the subject of future work. 
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APPENDIX A: CALCULATION OF r(a:) IN THE ANNEALED APPROXIMATION 

Let us tranform the partition function (|19p by exchanging the order of summation over confor- 
mations (uj) and disorder (\) and taking the explicit form of the sequence probability P(x), Eq ((4|) 
and later noting that the inner sum over all realisations of x can be written in the matrix form: 

(z N (f3\ x )) = £ P(x)£e xp • ^ £A + 0- - ^) eB ]f 

X uj ( i=l ) 



N 

^ ^ fxi e -|8Al '[ Xl£A+(1 ~ Xl)£Bl Y\Px i x- e" /3Ai ' [xi£A+(1 ~ Xl)£Bl 

u) x i—2 

Eft Ai r AinTW PAAW A l PABWg 1 \ I 1 
(fAW A \ fBWg 1 ) II £. A 4 • 1 



(Al) 



i=2 

N 



Pbaw a - pbbw b 
= ^(W Al f) T J](PW A ')-e 



i=2 



where we have defined wa = e 135 A , wb = e ^ £B , and 

'ft)' w =(7„° s ). .-(0- (A2 » 

The final result in Eq. (|A1|) does not depend on the sequence, i.e., we have effectively a ho- 
mopolymer partition function (but of a somewhat complicated form) . Since Aj can be either or 
1, each conformation will give us a product with the factors f T or Wf ^ (from the first monomer) 
and P and R = PW (from the monomers with i = 2,3...iV)ina certain order corresponding 
to this conformation. Unfortunately, the matrices P and R do not commute and we are not able 
to simplify this expression by, for example, rewriting the sum over w as a sum over the number of 
polymer contacts with the surface [l6l |. 

On the other hand, each adsorbed sequence (train), loop, or tail with the length n contributes 
the factor R", P n , or P", respectively, in the partition function (|Aip (if we do not consider 
the first monomer) and we can use the sequence GF method (or the grand canonical ensemble 
technique) suggested by Lifson [l8[ Introducing the GF of (Z n (B\x)) (or grand canonical p.f.) 
r(a0 = S^Li (Z n (f3\x)) an d applying Lifson's arguments [18| we obtain the result (fT!?j) for 
T(x). 



APPENDIX B: CALCULATION OF LOOP GENERATING FUNCTION 



In this Appendix, we briefly describe the calculation of the loop GF for the simple (6-choice) 
cubic lattice. First let us recall the "definition" of the loop. A loop of length n 
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FIG. 9: Decomposition of a loop on 6-choice simple cubic lattice. 



• starts and ends in the closest to the substrate lattice layer (site), 

• contains NO monomers in the surface layer, 

• the last, nth, monomer is followed by the step in the surface layer, 

• all the effects (statistical weights) at the beginning of a loop may be attributed to the pre- 
ceding adsorbed segment. 

First, let us notice that the minimum loop length for the given model is equal to 1. To calculate 
the GF of a loop 

oo 

T L {x) = Y J ^L{n)x n (Bl) 

n=l 

in the 6-choice simple cubic lattice, we have to take into account two possibilities (Fig. [5]): 

1. All the monomers of the loop are in the second layer (First "term" in the rhs of Figure [9]). 
The corresponding number of conformations is 

[z - 2) 71 - 1 

2. The first k monomers are in the second layer, then a step into the third layer follows. We 
decompose this contribution in three terms: The first k monomers, then loop starting from 
the second layer, which ends as soon as the first monomer reenters the second layer, and the 
rest. (Second "term" in the rhs of Figure [5]). The number of conformations is given by 

n-2 1=1 

J2 {z-2) k - 1 ili l {l)a L {n-k-l), n>3 

k=l n — k—1 

This gives us a recursion relation for f2i(n): 

n-2 1=1 

fl L {n) = (z - 2)"- 1 + Y (z~2) k - 1 n L (l)n L (n-k-l), n>3. (B2) 

k=l n — k—l 

By multiplying this equation by x n and summing from n = 3 to infinity, we obtain the following 
recursion relation for the GF Tl (x) (|Bljl 

r4«)-x-(,-2)^ = I ?^21 + T -^r| W . m 

Its solution is 

T L (x) = i- {l - (z -2)x - v/(l - zx)[l - {z - A)x]} . (B4) 
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